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COMMUTE  -  A  Computer 
Code  for  Noncommutative  Algebra 


1.  INTRODUCTION 

In  the  course  of  investigations  of  infrared  absorbing  frequencies  and  line 
strengths  of  triatomic  molecules  of  atmospheric  interest,  the  main  theoretical 
approach  has  invariably  involved  the  necessity  of  calculating  complex  commuta¬ 
tors  of  quantum  operators.  In  particular,  when  dealing  with  the  vibration- 
rotation  Hamiltonian,  one  encounters  commutators  and  anticommutators  of  angu¬ 
lar  momentum  operators  and  related  direction  cosine  operators  between  molecular- 
fixed  axes  and  space-fixed  axes.  This  report  presents  a  computerized  method  of 
generating  simplified  series  of  operators,  a  so-called  reduced  set,  from  given 
initial  commutators,  anticommutators,  or  general  operator  polynomials.  This 
computer  program,  which  implements  noncommutative  algebra,  is  called 
COMMUTE. 

2.  THEORY 

2.1  General  Formulation 

The  general  problem  is  one  in  which  there  is  a  commutator  consisting  of 
products  of  some  fundamental  operators,  whose  basic  commutation  relation  is 

(Received  for  publication  3  July  1979) 
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given.  Through  successive  application  of  this  fundamental  commutation  relation, 
the  original  commutator  can  be  simplified  into  a  series  of  anticommutators.  This 
can  be  expressed  symbolically  as 


C0[QJ,.Q'^1  *  CjlQ'j.Q'Jl  +C2{Q£.Q£}  +  ...  +Ci(Q!,Q,i'}  +  ...  ,  (1) 

where  the  O!  and  O'^  are,  in  general,  products  of  operators,  the  C.  are  either 
pure  real  or  pure  imaginary  numbers,  and  the  brackets  and  braces  are  the  stand¬ 
ard  notation  for  commutators  and  anticommutators,  that  is, 


l  A,  B]  =  Ab  -  BA 

and 

{A,  B}  =  AB  +  BA  .  (2) 


The  series  in  Eq.  (1)  in  practice  terminates  rapidly;  the  number  of  terms  is  pro¬ 
portional  to  the  degree  of  the  original  operators  OJj  and  Qq.  The  degree,  that  is, 
the  sum  of  the  powers  of  fundamental  operators  in  the  operator  polynomial  on 
the  right  side  of  Eq.  (1),  is  at  least  one  less  than  the  original  operator  polynomial 
[Oq,  Oq]  .  In  fact,  if  the  original  commutator  is  of  degree  n,  then  the  series  will 
contain  terms  of  degree  n-1,  n-3,  n-5,  etc.  (Of  course,  terms  may  be  missing 
altogether  if  their  coefficients,  CT,  are  identically  zero. )  Furthermore,  it  follows 
from  the  preservation  of  hermiticity  (or  antihermiticity)  that  if  CQ  is  an  imaginary 
number,  then  the  coefficients  of  the  anticommutators  are  all  pure  real.  Con¬ 
versely,  if  CQ  is  real,  the  coefficients  on  the  right  side  of  Eq.  (1)  are  pure 
imaginary. 

These  remarks  will  be  demonstrated  in  the  examples  that  follow. 


2.2  Diagonalizing  the  Hamiltonian 

The  problem  discussed  in  this  section  has  been  presented  by  Rothman  and 
Clough1  and  typifies  the  use  of  program  COMMUTE.  A  contact  transformation 


e1— '  H  e_i— ’ 


(3) 


is  applied  to  the  Hamiltonian  to  yield  a  transformed  Hamiltonian  expanded  in  a 
Taylor's  series 


1.  Rothman,  L.  S.  and  Clough,  S.A.  (1971)  J.  Chem.  Phys,  55:504. 
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00 


[S',  . . .  [S',  [S',  H]  1. . .  ] 
In) 


(4) 


H'(A)  =  Ji  +  Y, 
n=  1 


,n.n 
A  l 

n! 


where  X  is  an  expansion  parameter  and  S'  is  an  Hermitian  operator  (i  =  \Al). 
Utilizing  the  order  of  magnitude  arguments  given  in  Reference  1,  one  arrives  at 
expressions  for  the  "  rotational  Hamiltonian" ,  Eqs.  (31)  of  Reference  1.  S^'  is 
chosen  to  remove  certain  off-diagonal  terms  in  the  Hamiltonian  and,  as  such,  is 
itself  formed  of  angular  momentum  operators.  Thus,  expressions  for  the  rota¬ 
tional  Hamiltonian  contain  numerous  commutators  of  angular  momentum  operators. 
In  Appendix  B.  1,  two  representative  commutators  arising  from  the  contact  trans¬ 
formation  are  given. 

An  alternative  scheme  for  calculating  the  rotational  commutators  has  been 

2 

presented  in  the  work  of  Amat  et  al.  In  their  approach,  the  contribution  of  a 
chosen  operator  term  is  given  as  a  function  of  commutators  of  various  directions. 
As  such,  their  procedure  is  the  inverse  of  the  one  discussed  here. 

2.3  Reducing  Operations 

The  rotational  Hamiltonian  through  fourth  order  of  a  planar  molecule  in  the 
XZ  plane  is  of  the  form 


H™  =  h  +  AP2  +  BP2  +  CP2  +  D  (P  ,  P  }  +  Y 
—  v  — z  — x  — y  1 — x  — z  ’ 

afad 


p'  c  P  P„P  P, 
ta/3->6  —  o-p — : — 6 


ofiyit  n 


P  P,jP  P.P  P 
op-yScn  —a—fi—y— 6—e—r) 


(5) 


where  hy  represents  the  vibrational  contribution  to  the  Hamiltonian;  A,  B,  C,  and 

D  are  coefficients  of  the  quadratic  terms  in  the  angular  momentum;  the 

are  coefficients  of  the  quartic  terms  in  angular  momentum;  and  the  $'  -  .  are 

B  affydcri 

coefficients  of  the  sextic  terms  in  angular  momentum.  For  the  case  of  ortho¬ 
rhombic  symmetry,  thereare  15  nonvanishing  p  ^  and  105  -A  great 

amount  of  the  simplification  can  be  obtained  in  the  Hamiltonian  by  application  of 
the  commutation  relations  of  angular  momentum  to  these  coefficients.  The 


2.  Amat,  G.,  Goldsmith,  M. ,  and  Nielsen,  H.  H.  (1957)  J.  Chem.  Phys.  27:838. 
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resulting  reduced  form  of  H^7  for  asymmetric  rotors  of  orthorhombic  symmetry 

—  3 

thus  contains  only  19  coefficients  (Kneizys  et  al  ): 


H 


tt  =  h  +  A  P2  +  B  P2  +  C  P2 
v  — z  —x  —  y 


4  4  4 

+  t1p^  +  t2pJ  +  t3p4z 
+  T4  [p2  P2}  +  T5  {P2,  P2}  +  T6(p2,  p2) 


+ 


.  P  +  *„  P  +  *,  P  +•  *,  {P  ,  P  I 
l— x  2— y  3— z  4'— x'— y’ 


+  *5  (Py-4t  +  «4  p4z>  +  *7  &  4} 

+  {p2,  P4f  +  4>„  (P2,  p4}  +  $.r.(P2  P2  P2  +  P2  P2  P2) 

8  1 — z  — x  '  9  1 — x* — z  ’  10  — x — y—  z  —  z — y  — x 


(6) 


2.4  Minimal  Set  of  ()|M*rators 

4 

It  has  been  shown  by  Watson  that  the  coefficients  of  the  Hamiltonian  of  Eq. 
(6)  are  not  all  linearly  independent  for  an  asymmetric -top  molecule.  This  gives 
rise  to  highly  correlated  coefficients  when  an  attempt  is  made  to  estimate  the 
coefficients  in  conjunction  with  a  least  squares  approach  with  the  experimental 
data  (Flaud  and  Camy-Peyret,  Rothman  and  Clough  ).  The  indeterminancy  can 
be  further  removed  by  successive  unitary  transformations  on  of  the  form 

H  =  e-H^  e‘1^  .  (7) 

Successive  transformations  of  this  type  again  lead  to  a  series  of  commutators 
analogous  to  Eq.  (4).  The  choice  of  elements  where 


1  1 


s  (PpPqP'  +  P‘  P4?^) 
pqr  -x— y-z  -z-y-x 


3*1  t>P\ 


pqr 


(8) 


3.  Kneizys,  F.X.,  Freedman,  J.N.,  and  Clough,  S.A.  (1966)  J.  Chem.  Phys. 

44;2552. 

4.  Watson,  J.K.G.  (1967)  J.  Chem.  Phys.  46:1935;  (1968)  J.  Chem.  Phys. 

iS:4517‘  — 

5.  Flaud,  J. -M.  and  Camy-Peyret,  C.  (1973)  Mol.  Phys.  26:811. 

6.  Rothman,  L.S.  and  Clough,  S.A.  (1975)  Determination  of  Valence  Force 

Constants  for  Water  from  Vibrational  Data.  Paper  Presented  at  the 
Thirtieth  Symposium  on  Molecular  Structure  and  Spectroscopy,  Ohio  State 
University,  Columbus,  OH. 
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is  arbitrary  and  can  be  made  to  reduce  the  Hamiltonian  to  contain  only  (n  +  1) 

independent  terms  of  total  degree  n  in  the  components  of  total  angular  momentum 

for  each  even  value  of  n.  The  standard  choice  is  to  remove  the  coefficients  of 

operators  such  that  the  matrix  elements  in  the  symmetric-top  basis  satisfy  the 

selection  rule  AK  =  0,  ±2,  that  is,  the  reduced  Hamiltonian  matrix  in  a  given 

7  8 

vibrational  state  is  at  worst  tridiagonal  (Yallabandi  and  Parker,  Benedict  et  al  ). 
This  leads  to  the  Watson  reduced  Hamiltonian 


»w  =  h000  +  h100  P2  +  h010 4  +  h001  <4  -  4] 


+  h200^  +  hao4  Z\  +  h020P4z  +  h101  |P,2(p2  -  p2» 


,2,„2 


+  h011  <4  <4  -  4)}  +  h300— B  +  h21044 


+  h12o44z  +  h030P«  +  h201  (P4(p2  -  P2)} 


^  hul  {P2  P2Z.  <P2  -  P2)}  +  h021  {P4.  (P2  -  P2)} 


(9) 


Eq.  (9)  is  given  through  fourth  order,  although  many  terms  of  high  order  in  angu¬ 
lar  momentum  can  be  preserved  in  the  transformation  Eq.  (7)  so  that  the  equality 
is  preserved.  Notice  that  the  choice  of  the  parameters  Spqr  has  been  made  to 

9  P 


2  2  r  n" 
eliminate  any  powers  of  the  operator  (P^  -  P  )  greater  than  unity. 


2.5  Direction  Cosine  Operators 

Another  example  of  the  use  of  the  commutator  reduction  is  seen  in  the  calcu- 

9  10 

lation  of  the  dipole  moment  expansion  (Clough  et  al,  Flaud  et  al  ).  This  prob¬ 
lem  requires  the  calculation  of  commutators  of  the  type 


K?£p?pz  +  £44.  .  <9 io> 

where  the  indices  p,  q,  r  are  powers  of  angular  momentum,  and  9^  is  the  direction 
cosine  between  the  o-molecular- fixed  axis  and  the  space-fixed  axis.  The  com¬ 
puter  code  COMMUTE  generates  series  of  reduced  operators  for  the  general 


7.  Yallabandi,  K.Y.  and  Parker,  P.M.  (1968)  J.  Chem.  Phys.  49:410. 

8.  Benedict,  W.S.,  Clough,  S.A.,  Frenkel,  L. ,  and  Sullivan,  T.E.  (1970) 

J.  Chem.  Phys.  53 :2565. 

9.  Clough.  S.A.,  Beers,  Y. ,  Klein,  G.P.,  and  Rothman,  L.S.  (1973)  J.  Chem. 

Phys.  59:2254. 

10.  Flaud,  J.-M.,  Camy-Peyret,  C.,  Mandin,  J. -Y. ,  and  Guelachvili,  G.  (1977) 
Mol.  Phys.  34:413. 
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term  (10)  in  a  fashion  similar  to  the  results  from  the  previous  examples.  The 
most  general  form  of  commutation  operations  involving  the  direction  cosine 
operators  occurs  when  performing  a  Watson-type  transformation  on  the  dipole 
moment  expansion.  The  form  of  the  commutator  is  then 


[(PpPqPr  +  Pr  PqPp) 
— x — y _ z  — z — y — x 


(p8p*pu*  +®  p'Vp8)] 

— x  — y  — z  -a  —a—  z— y— x 


Examples  of  the  application  of  COMMUTE  to  terms  (10)  and  (11)  are  given  in 
Appendix  B.  4  and  B.  5. 
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Appendix  A 


The  Code  COMMUTE 


Provided  in  this  Appendix  is  a  complete  listing  of  the  program  COMMUTE. 
Fundamental  commutation  relations  are  built  in  for  the  angular  momentum  opera¬ 
tors  and  direction  cosine  operators.  To  distinguish  the  two  types  of  operators,  a 
general  6-vector  was  chosen  with  the  correspondence:  U  —  ®  ,  V  -*  ®  ,  W  -*  ® 


general  6-vector  was  chosen  with  the  correspondence:  U  —  V  —  ®  W  — 

X  -  P  ,  Y  -  P  ,  and  Z  —  P  . 
x’  y'  z 

The  first  input  card  indicates  the  number  of  terms  in  the  initial  polynomial 

that  is  to  be  reduced.  However,  the  program  automatically  takes  each  initial 

term  and  forms  an  anticyclic  pair  to  it.  Thus,  for  the  commutator 

[P  ,  P  ]  =  P  P  -  P  P  ,  it  is  only  necessary  to  consider  one  term  P  P  .  The 
xyxyyx  x  y 

operator  polynomial  [ { P  ,  P  },  P  ]  should  be  considered  as  two  terms  for  the 
x  z  y 

program,  P  P  P  and  P  P  P ;  anticommutators  should  be  expanded.  This  fea- 
x  z  y  z  x  y 

ture  of  doubling  terms  was  undertaken  since  the  primary  use  of  the  program  re¬ 
quires  commutator  or  anticommutator  reduction.  Single  terms  can  nevertheless 
be  handled  as  can  be  seen  from  the  example  in  Appendix  B.  There  is  also  the 
option  on  the  first  input  card  to  allow  the  program  to  drive  the  reduction  in  a 
cyclic,  that  is,  UVWXYZ,  or  anticyclic  direction  (ZYXWVU),  thereby  facilitating 
tests  for  axis  invariance. 

The  second  group  of  input  cards  (equal  to  the  number  of  initial  terms)  indi¬ 
cates  the  number  of  operators,  the  coefficients,  and  the  directions  and  exponents 
of  the  operators.  The  sample  data  deck  accompanying  the  examples  in  Appendix  B 
should  suffice  to  explain  the  method  of  reading  in  the  terms. 
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The  program  works  on  two  "levels"  at  any  time,  a  level  being  all  terms  of  a 
given  degree  in  exponent  sum.  The  operators  are  scanned  for  their  directions. 

If  the  directions  are  different,  then  the  commutation  relation  is  applied,  if  they 
are  the  same,  then  they  are  multiplied  (exponents  added).  The  action  of  the  com¬ 
mutation  relation  creates  a  term  at  the  next  level.  Terms  at  the  second  level 
are  scanned  for  identical  terms  and  cancelled  when  the  coefficients  are  zero. 
When  a  level  has  been  completely  scanned  and  condensed,  it  is  written  as  anti- 
commutators  if  the  original  term  or  terms  was  either  hermitian  or  antihermitian. 
The  process  thus  is  one  in  which  commutators  are  always  "driven". 
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«  ART  TCOPMUTI  TOF  S)  THROUGH  THE  US*  OF  COMMUTATION  RFLATIOKS. 

FOR  Tvip  ARRAY*-  JH  PLANY  COMPCNt 

F  T°rT  IPOtF  Or'K-'FTES  LEVEL,  S.CONDUERH,  TPJRC  IPOS  IT  1  C  K  CF  CFERATCI 
IFO«HA»FlfV,I  TEFY)  S  NOHREP  OF  OP;PAT0RS  IK  ITE  PH. 

TTEFt-LELFVi  ■=  TrTfl  NVNFER  OF  TFR1S  IN  LTV. 

rcQCFJM  Mil  L  Trt,i(i(ic  IF  jj  PNCCUNTFRSt  (1)  TMFLVS  OR  HCPE  CFERATOR: 
H  6  TCC-H,  OR  I?)  ?nO  OB  “ORE  TT->MS  AT  A  PAOTICULAR  LEVEL. 

COHSON  T  *  T  e  I  t  ,1G9,11)  ,  1  E  XP  ( 2 , 1 R  9  *  11), COcFf  (7,199) ,  I F CSMA  X ( 2, 1 99) 
CIHFKS1  )H  ITCF  cu?)  ,  PIR(6>,  X(11>,  X2(ll),  EPS(6,F>,  ICTEL(F,6) 
lHTYPdi),  )Fir(P),  tETFEdl),  F  FS  VEC  ( 36 )  ,  ICVFS  (36) 

CATA  <riPn),Trl,F,)/lHH,lHV,luW,lHX,lHY,lPZ/ 

CATS  (HOIG(l)  ,I*l,9)/lHltlH?,lHT,lH4,lH9,lHf,lHI,lHP,lHS/ 
rire  mp.mlf ,mi  f  ,hpl  ,hhi, hr/ihf,  ih  «,  in) ,  ih«,  ih-,  ih  / 

rOHH(iIATI0V  B  rl  FT  7  P*'S  I 
OAT  A  FF  X  V  ~  C  / 
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0  9 
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9 

1  9 

r  9 

c 

U  9 
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COP 

1C 

COH 

11 

COP 

12 

COH 

13 

COH 

14 

COM 

15 

COH 

16 

COH 

20 

COH 

21 

COP 

22 

COP 

25 

COP 

30 

•  COP 

40 

COP 

50 

COM 

65 

COP 

70 

COH 

so 

CCH 

9C 

COM 

100 

CCP 

no 

COM 

120 

COH 

130 

CCH 

140 

cop 

1  50 

COM 

160 

COM 

170 

COP 

ltli. 

no  IP 

co  n  jd.f 

INP  -  J*4  *  (  1-1  > 

FFS  (  I,  J)  *T  F  svt  C(  INC) 
TOT  CL  (I  ,J>  cICVEC  <nc> 
CCK'TIHJE 


10 
c 

C  (1)  F  "AC  IM  V'i»prF  OF  INPUT  TERM 

c 
c 
c 
c 


CCP 

COP 

COH 

CCH 


(17)  ARC  HHfTMEF  T frhs  ARE  TC  9E 
ORTVES  CYC'.Ir*LLY  OF  AM1ICYCLIC1LI  Y  (F&.l); 

(7)  P-.Af  NU-'o-F  o'  CPfOATOPS  (17),  CO'FFICIENTS  (7FF.1),  ANC  CIFECTICNCOH 
ANO  CYPONFPT  OF  rrrv.  pocoator  —  FOPHAT  1 1 2,  ?F  6 . 1 , 1 1  ( 1 X  ,  A 1 , 1 7 ,2X  )  )  I  COP 

COM 


310 

32. 

3X0 

TAG 

35P 

36  0 

j  7 ; 


L "  V-  1 

COM 

1-ViT z 

COM 

390 

20 

e'lr  900,  ml  (1)  ,CYCl.  5 

COM 

kOO 

rrrsiPM";p"i  i" 

COM 

kio 

t  f  ( r  vf  l  r  .if.  -j.n  cvc«i  «mi\nticyclip 

COM 

k  JO 

vr¥  i  >*?•  rf  *»«u  n 

COM 

k  10 

ML*  I  Tf  *  1  ) 

CCM 

k  k  0 

CALI  snvir  (MMf) 

COM 

k5P 

p^TT'T  ,  TT^r  f  ^vr 

COM 

k  60 

TF  lITl.^.n  c  T  OP 

COM 

k  70 

co  f o  r  tt«  ,*  u  , ? 

COM 

fcJC 

*rAr  <r*  ,  *A*#CPFFF  (1,1T>  #COcFF  n 

,1  Ttl  1  ,  CM  IP) ,  IFYF  (1  ,IT  ,IF)  ,  IP=  1  ,cc» 

k  90 

1HAV  > 

COM 

500 

TcPSM*X<t  f  yn  *  I  POSH  AX  (  1  ,  IT  4  1  )*MAX 

COM 

51. 

ro  ’0  I  * 4  * M  f  x 

COM 

6  2  0 

ro  »fl  J:1,f 

COM 

53: 

TF  (PTP(  J>  .eo.x  CD)  I  01  ^  ( 1  *T  T  #1  )  * 

J 

COM 

5k. 

xr 

row^TNUc 

COM 

660 

r 

COM 

660 

r 

A MT  CMAT  IC  Al  L  *  cv*1FT*IZ£  INP(Jt  Tr*HSI 

COM 

670 

r 

COM 

55  0 

T  «ll*»  p 

COM 

690 

ro  41  |*1  , 

COM 

6  j  . 

TSUm*  ycim^Tr  yr  ( i  ,  IT  fT) 

COM 

610 

J*Mf  y*i  -T 

COM 

620 

niMi. «ni®  citiT 1 1) 

COM 

630 

Trxrci#tT*1 ,T»  *7px" (lt  tt , J> 

COM 

6k" 

M*7FX?M  .17, M 

CCM 

650 

K  j*  yr yp(  i  ,TTM  ,  I» 

COM 

660 

4  =  XC  (II  *4*71"  (NT> 

CCM 

670 

wpyp’ a  \  *  «r jc<  nji 

COM 

680 

40 

X’( T 1 *X ( J1 

COM 

690 

r 

COM 

7.  . 

c 

TF.T  fnc  -1  FOR  kO>l-M£PMlTl»K  Iff UT  CP6ft»TCk,  COM 

710 

c 

fl  f  np  urijMlJHl  Of  f  RflT"F,  *1  FOk  »MTI 

-HF^HITIAN 

OPt  PA! OF  1  COM 

720 

c 

Com 

731 

I"#  I  sff-'f  (J,lT«l)/fOeFF  (1,IT)*1, 

1 

COM 

7k0 

ic« rs Tr'i“-?»  (tsiim/’) 

COM 

750 

IMF  F  Ms  1 

Com 

760 

IF  ( (!P»»  .  FC  .  M  .AND  .  CI»«T.FO.  «  )  ) 

lHtPM  =  0 

COM 

770 

I f  i(TF«p.Fc.ri,**io,(iP»T.Fij.?n 

IHFPKrO 

COM 

7BC 

IF  (T9Ar,“",l)  IMFOMs-l 

COM 

79r 

r 

COM 

moo 

f 

PRIM  ISFirT  T 

COM 

6  10 

r 

COM 

820 

tf  i (m*»,  ro, i)  .or.  i  lMrPM.io.  -  in 

r,0  TC  *0 

COM 

0  3i 

PJT»T  Ok",  (M|  V  f  (  IP  )  ,««  ,  10:  1  ,  M«  Y) 

♦  H? » 

(HFXF2I  IF  )  ,Hfi  ,  IFM  ,  “l>X)  COM 

6  kO 

COM 

0  50 

IF  UM»T."Q.’>  vnirf'sHTL 

CO“ 

06. 

r«Li  mr’  <cffff  <i  ,iti  ,tc,no> 

COM 

670 

tf  iMr.to  n  ffi«it  ‘•to,  ic,Ma,HiP 

,HO,  (HP  ,^o 

,IP*1  ,M«XI,HB,XSJf.N,  COM 

6  n  0 

1  CM"  ,HP,  T"s  1  ,90  Y  )  ,MI)  ,Ml>r 

Com 

0  90 

IF  Al.T.  c  N0.Lp.F*4>)  PCIKT  9  ?P,  1C 

,N0,HB,HIF  ,HB  ,  (MP ,HB,  COM 

9"0 

1TP*  1  . »  °tvr  IGN,  Tr*1  »HAX) 

COM 

610 

TF  PkTNT  93d  CO?ff(i, 

it  )  ,ne,Hi  p 

,  HP ,  tHP,M6,IF>l,M«X)  ,MBC0M 

Q  20 

«  f  »SIRF.  f  TF*1  f  *AX) 

COM 

930 

P3I»T  $uP,  (V r TF  )  ,HFt  IPs  1 # MAX)  fMP 

,MP,MO,  (V’  (IF)  ,Ma,  Ifrl,M/>>)  COM 

9kt 

18 


50 


GC  10  6' 

PRINT  940,  (MrXP(IP) , H  B , I®=  1 ,  MA X) 

C=COEFF(1,TTMCOEFF(l,IT*l> 

call  f»*:t  (C.ic.nd) 

IF  (N0.cn. r>  PRINT  910,  IC,  H9 ,  H0  ,  MB,  (  HF, HP ,  IP*1  ,MAX  I 
IF  (In6.GT.CI  .ANO.(NO.LE.64)  >  POINT  920,  IC,NQ,HB,H8,FP,  (HP, HP, 

1  IP=1  ,MAX> 

IF  (NO.Ji.'l.)  POINT  970,  C,MP,HB,HB,(HP,HP,IP«1,N»X) 

POINT  940,  (X(IP)  ,MP,IP*1,MAX) 

60  CONTINUE 

POINT  90** 

C 

C  INOFX  LFVELS,  ti-rrs,  AMO  TEST  POSITION  OF  OFEPA  TORS  (SHIFT  LEVELS  TO 
C  SAVE  FTOPACd  I 
C 

L  =  0 

70  CUM 

if  (l .n.i)  gf  to  oo 

ITL*ITfOML(LEV2> 

IF  (ITL.FO.f)  G  C  TO  70 
ITFPHL'Lf  V)  =  TTL 
00  FS  IFs1,TtI 

coeff(i=:v,tt*scpfff  (lev?, id 

max=IP09«»7(LF  V7,  IT) 

IPOSMAX  (Lcv,  IT  )  =  «AT 
PO  60  rP*1,MAX 

IOIP(Lf  V,  IT,  IF  )  =  IOIP(  LFV2.IT,  IP) 

0  0  TEX P  (LEV  ,  TT,  IP!  =  IEXFILFV2.it,  IP) 

90  ITFONL(LE'<7)*ITfOM*ITEPH2sn 

100  ITE  PM=I Tr  ®F  ♦  1 

IF  (ITfOML(L'V).EQ.ri  GO  TO  70 
IF  (IT?  PN.FT.ITFRMMLEV)  )  GO  TO  740 
110  IFOF*0 

120  IPOS*inF*1 

IF  (IFOS.GF  .  IFP?:nax(LFV,  ITFR1)  )  GO  TO  7qo 
N=TPOF*l 

if  (101°. (LFV.ITEOH.N)  ,EO.InI9(LEV,ITEPN,IPOS>)  GO  TC  130 
Ie  (CVCL F  ,LT.  -0.1)  GO  TO  171 

IF  <ICd(l  FV,  I  TERN,  IPOS)  .GT.  I0IR(LEV,  ITEPM.N)  )  GO  TC  160 
GO  TO  »7 P 

121  TF  (ICTP(LrV,  ITFOM,  IPOSI.LT.  IPIR(LEV,  ITEOH.Nl)  GO  TO  160 
GO  T0  (20 

C 

r  CONTRACT  APJ  ■  nf  nt  OPER« TO»E  HAVING  IOENTICAL  DIRECTIONS,  IN  LEVA 
C 

13"  1FXP(Lf  V,  TTrop,  IPOS)*IEXP(LF'/,ITFRN,IPOS)  *1  EXP < L EV , ITERM  ,N  ) 

NAXrlPCSNAX (L>  v  ,  IT? FN)-1 
IF  (MAT.'.t.N)  GO  TO  150 
ro  540  T*N,  MAX 

ICIP(U  V,  TTF  PM,  I  )  =  TriR  (LEV,  I  TF  RN,  I»l) 

140  iFXPdfv,  nrRN,i)=ifXP(LEv,iT£RM,i*i) 

150  IFOSNAr  (LFV,ITFON)  =  pAX 

CC  TO  110 

c 

C  CONFUTE  TWJ  OPFOAIORS  IN  LEVI 

c 

160  I  A*  IDIR(lrV,ITFPN,TFOS) 


CON  950 
COM  960 
COM  970 
CON  960 
CON  990 
COM  1000 
CON  1010 
CON  1020 
COM  1g  30 
COM  1040 
COM  1050 
CON  1060 
CON  1070 
CON  1.60 
CON  1090 
COM  1100 
CON  1110 
CON  1120 
CON  1131 
CON  1140 
COM  1150 
COM  1160 
CON  1170 
CON  1160 
COM  1190 
COM  1200 
CON  1210 
CGN  1220 
COM  1230 
CON  1240 
COM  1250 
COM  1260 
CON  1270 
COM  1260 
COM  12  90 
COM  1300 
CON  1310 
COM  1320 
CON  1330 
COM  1340 
COM  1350 
COM  136t 
COM  1370 
CON  1760 
COM  139C 
CON  1400 
COM  1410 
COM  1420 
CON  1470 
CON  1440 
CON  1450 
CON  1460 
CON  1470 
CON  1460 
CON  1590 
CON  1500 
COM  1510 
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t3sI0TMlr'', ITERH, N) 

IEXF(IEV,TTE9M,  JPOS>*IE<B(LFV,  ITERH, I  FOE )-l 
IpORmAMi  EV,I’  E  RM)*1P0SHAX(17V»ITERN|  ♦» 
hAXsIPOSHAt  (Lf  V  ,ITFPM> 

IF  (PM. IF. 11)  rc  Tr  170 
PRINT  ,  MAy  , l, TERM 
CO  TO  ?0 
170  LIMsMAX-" 

ro  loo  t*i ,i in 

Jt>*(X*l  -T 

iriP(LEtf,TTE01  ,  J)*iriR(LFV,rTERN,J-?l 
16  0  T'XPd  E'MTrR>*.  J)  *IFXP<LEV,IT?»Mf  >71 

IEXP(UV,  ITFRN,N*2>rIFXP(LrV  ,1  TERM,N)  -1 
IETP  (l  CV,  TTFRH,n*1>  *IF  XP  ( lc  V,  I TFRM,  N)  *1 
C 

C  ELIMINATE  7E®''-EXPPNFNTR  IN  LeVI 

C 

N»»Nt2 

KCAMsN 

TF  (IFX',n.n/,lTEBM,N2).GT.  D  00  TO  21 T 
MAX*  IPO'P'X  (Lr  v  ,  ITFPM)  -1. 
if  (Hex.  lt.n>)  oo  Tr  ?oo 
rr  190  t=M?,ifx 

inikdtv,  TTFRH ,  I)  «I PIP (LEV,I TERN, IH) 

190  IEXP(lFV,IT'RM,I>*tFXP(lEV,IT£RM,T*l» 

200  T  POSH  AX  (l  rV,  ITERH)  =  MAX 

210  IF  (IE*P(LF  v,ITFRM,IPOS)  ,GT.  0)  GO  TO  270 
KOAHs  NG»"-1 

HAXsIPOSHAXCLrV,  ITERM1-1 
90  220  1  =  TrOS*  M6  X 

IPIPdf  V,ITFRH,I)*iriR(LEV,ITERH,I»l> 

220  IEXFUEV  ,ITERK,I)=IEXP<LEV,ITERH,I*l) 

IROSMAX  (  E  r  V,  IT E  PH)  *  hA X 
230  IF  (FP>  (r  A, IP)  .EO.O.)  ro  TO  110 

r 

r  CREATE  TERM?  IN  LE  V2  <s  LtV*l)l 
C 

1TERH2:  I»^FH?*i 
IF  <ITf«>H?,ic.  jgq)  GO  TO  »40 
PRINT  RAO,  ITFRM2,l,IT£PM 

gc  to  n 

240  ITFRHLtLrV?)sITFRH2 

I0OSMAX(LCV?,TTERM?)*IPCSNAX  (l*V,  ITERH) -1 
TRsN'GAH-1 

if  (1F.LT.1I  ro  TO  260 
no  ?ro  t=< , ip 

iriPCLf  V2,IT  =  FH2,I)sIOIF(LfV,ITERH,I) 

260  lEXPdFV’.ITFf  H2,I)*IEXP(LCV,  ITERH,  II 
260  rcFFFUEV»,ITFPH?)*tPF  ( I A , 10 > *COEFF  <L  EV  ,  ITERM) 

I0IP(LEV»,ITEFH2,NGAH)*ICTBLCIA,I9I 
IFXP(lfV2,TTcrH2tNr.FH)*l 
HAXsIPOSi'X  <IEV2,ITFRH2> 

N1  *  N  C  AM  ♦  1 

TF  (MAT.  IT.  Hi)  r.O  TO  2  89 

ro  27C  Tsmj,m«x 

I0IF(«V»,ITFFH?,I)  sIOIP  (LEV,  ITERH,  I»l) 

270  IEXF(IFV7,ITEFH2,T)  =  IFXF(LTV,ITERM,  1*1* 


COM  1520 
CCiH  1530 
COM  1540 
COM  1550 
COH  1560 
COM  1570 
COM  1500 
COM  1590 
COM  1600 
COM  1610 
COH  16  20 
COH  1630 
COH  1640 
COM  1650 
COM  1660 
COH  1670 
COM  1660 
COM  1690 
COH  1700 
COH  1710 
COH  1720 
COM  1730 
COH  174C 
COM  17  50 
COM  1760 
COM  1770 
COH  1700 
COH  1790 
COH  1600 
COM  1610 
COM  1820 
COH  1630 
COM  1840 
COM  1650 
COM  I860 
COH  187t 
COH  I860 
COM  1690 
COM  1900 
COM  1910 
COH  1920 
COM  1930 
COM  1940 
COM  1950 
COH  1960 
COM  1970 
COM  I960 
COM  199C 
COH  2000 
COM  2010 
COM  2020 
COM  2030 
COH  2040 
COM  2050 
COM  2060 
COH  2070 
COH  2060 


c 

COH 

2090 

C  CONTRACT  A0 JAG(  NT  OPERATORS  HAVING  IO£NTIC*l  DIRECTIONS  IR  L EV2I 

COM 

2100 

C 

COH 

2110 

260 

I*  (HAT.  IT.  21  ro  TO  310 

COH 

2120 

l»AX*HAy 

COH 

2130 

LI*  1 

COH 

2140 

ro  300  T*»,HAX 

COH 

2190 

11*11*1 

COP 

2160 

IE  (LHAT.LT.LT)  GO  Tq  J10 

COH 

2170 

IE  (IOIXLfV?,ITFRH7,Ll).NE.IDIPClEV2,  ITERH2.LI- 

1))  GC  TC  300 

COH 

2100 

If  XPIlf  V’,  TTEF  H2,  LI  -1)  *IE*P(LE  V2,ITERH2,LI 

-1»*IEXP(LEV2,ITERH2, 

L  I)  COH 

2190 

LPAX*LPA»-1 

COH 

2200 

IPOSHAI  (irv?.ITERH2)*LHAX 

COH 

2210 

IE  (LHAX.LT. LI)  GO  TO  310 

COP 

2220 

ro  J90  1*11, LPAX 

COH 

2230 

10  IP(  Lf '/?,  TTCR  H2,J)-I0IR(LEV2,ITERH2,J*1) 

COH 

2240 

290 

IfT  f  ( LI  <•?,  IT*  T  H2,  J)*  IE  XP  (LEV?,  ITERN2,  J*l) 

COH 

2250 

LI*  LI-I 

COH 

2260 

TOC 

PONT  1MJC 

COH 

2270 

C 

COH 

2200 

C  AOO 

IDENTICAL  TERHE  IN  LEV2I 

COH 

2290 

c 

COH 

2300 

J10 

Tf  II  TTD-?.L  E.  1)  GO  TO  110 

COP 

2310 

PAX*  IPCSHA  y  (LFV2,  IT ERH?) 

COH 

2320 

»*ITEOH’-l 

COH 

2330 

no  T 3 9  I*', < 

COH 

2340 

If  IIP03H*XILtv2,I)  .NF."AX)  GO  To  330 

COH 

2350 

no  T20  J*1 ,«AX 

COH 

2360 

IE  (ini’ll  fv’,1,  J).NE.iriR(LEV2,ITERH2,Jl) 

GO  TO 

330 

COH 

2370 

If  (TE*»(l  CV?,I  ,  J).  HE.  IfX°(trV?,ITERH?,JI) 

GO  TO 

33  C 

COH 

2360 

520 

CONTINUE 

COH 

2390 

CCEFF  ILEV2,T)*C  OEFF  (l  F  V  2  ,  T  )  *COEF  F  (l  EV  2  ,  IT  FRP  2) 

COP 

24  00 

l’E  FH?» IT  E  LV2- 1 

COP 

2410 

ITEPNLILEV?)*  ITERN2 

COH 

2420 

CALL  CA VrEl  II  *V2,ITERH2,1,1TERHL(LEV2>) 

COH 

2430 

GC  to  1 10 

COP 

2440 

3  50 

CONTINUE 

COP 

2450 

GO  TO  1  10 

COH 

2460 

c 

COH 

2470 

C  TES 

T  EOF  AN!  T-rOHFUTATORS  OR  GOHTUT  A  TOR  S  t 

COH 

2400 

c 

COH 

2490 

»4fl 

TF  (IHCRH.rc.-l)  GO  TO  460 

COH 

2500 

1 SDH*  0 

COH 

2510 

PAX*  IP’GM»  X  (L  f  V  ,  1  ) 

COP 

2520 

no  'so  t * ’ , ha y 

COH 

2530 

350 

TSUH*  IS  'I**  JEyf  (LEV, 1,1) 

COH 

2540 

IOAP«I3'l*-?*(T*HH/2) 

COH 

2550 

IF  (IlfoH.rO.lPAR)  GO  TO  4f« 

COH 

2560 

1tL*  I  Tf  R*M  1 1  -  >  » 

COP 

2570 

If  (ITL.lT. IIP)  GO  TO  360 

COH 

2560 

°RI t  T  90i-  ,  in  ,  L  ,  IT  (  RH 

COH 

2590 

GO  TO  ’0 

COP 

2600 

JfO 

OO  tji)  tt  =  i  ,  I  t  l 

COP 

2610 

PAX*IP0?H«X(L;  V,IT) 

COP 

2620 

l°OSHAX  II  *V,  11  ♦  IT  L)  *NAX 

COH 

2630 

err FF  (LEV,  TT)r  £  .S*rrEFE  (L- V,  IT) 

COP 

2640 

CnEFF  lit  V  ,n*lTL>  =-ro:FF  (LEV,  IT) 

COH 

2650 

*  ™  ^  -*  ■ - UIIJ,™  - - - -  -ll  . .  . . — 


M  370  Trt,1»7 

COM 

2660 

J=HAX  +  l-I 

COM 

2670 

IEXf  (UV,TT*ITl  ,J)=IEXP<L€V,  IT,  II 

COM 

2680 

370 

laiPILEV,  H  +  IU  ,  J):  IOIR(LEV,TT  ,1) 

COM 

2691 

380 

CCNTINVE 

COM 

2700 

TIE  PMsF  7|  +1 

COM 

2710 

ITFPMLCLiV)=2»ITL 

COM 

2720 

go  to  lir 

COM 

2730 

C 

COM 

2740 

C  *00 

IDENTICAL  lEPMf  IN  LEVI 

COM 

2750 

c 

COM 

2760 

T9  0 

IE  (ITE^i.FO.l)  GO  TO  440 

COM 

2770 

max*ifcsmiy<lfv,  itepmi 

COM 

2780 

K=ITERN-l 

COM 

2790 

PO  430  It  1 , K 

COM 

2800 

IF  (IPOENAXCLFV.II.NE.MAX)  GO  TO  430 

COM 

2610 

PO  400  Jtl  ,«AX 

COM 

2820 

IF  ( IPTRJ  lf  V,  J  ,  J>  .NE.IDIR(LEV,ITERN,J>  >  GO  TO  430 

COM 

2830 

IF  (IFX  =  «LFV,I,JI.NF.IFXF<LFV,ITE3M,J)>  GO  TO  430 

COM 

2840 

400 

CONTINUE 

COM 

28  5G 

CCEFFdEV  ,H*roEFF(LEV  ,I>  ♦COEFF  (LEV.ITFBM) 

COM 

2860 

ITFRML(LE'M  =  ITFRnLUEV> -1 

COM 

2870 

ITL= ITE?mL (LEV) 

COM 

2880 

IF  (ITL.Lt.ITE PM)  GO  TO  4?0 

COM 

2890 

CO  410  IT  =  1TERM,ITL 

COM 

2900 

LIMtIPOSMA>  (LFV,IT»1> 

COM 

2910 

C0F  FFILFV,IT)  =COEFF  (LEV,IT  +  1> 

COM 

2920 

I°OEMAI  (LEV ,  IT )  =LIM 

COM 

29  30 

PC  41 C  I°=  3 , L I M 

COM 

2940 

IOIP<LEV,TT,irj=lPlP(LEV,IT  +  l,IPI 

COM 

2950 

410 

IEXF  <LE«,  T7,  ID*IEXFtLFV,IT  +  tf  IP1 

COM 

2960 

42C 

ttfpmsitepm-i 

COM 

2970 

CC  TO  43 " 

COM 

29  8  0 

430 

CONTINUE 

COM 

2990 

CC  TO  10" 

COM 

30  00 

C 

COM 

3010 

r  CANCEL  7ER0-CO  EEpT  Cl  ENT  TERMS  IN  LEVI 

COM 

3020 

c 

COM 

3030 

uun 

1  =  1 

COM 

3040 

450 

CALL  CA‘|CM  (1  FV,ITERM,I,mRML(LEV>> 

COM 

3C50 

GC  TO  100 

COM 

3,60 

c 

COM 

3070 

C  VDTTF  AM 

COM 

3080 

r 

COM 

3090 

460 

ITL=ITf5ML  n<ru) 

COM 

3100 

lc  (ITl.-c.n  CO  TO  70 

COM 

3110 

I'  (L.ET.ll  G?  TO  470 

COM 

3120 

LL=L-1 

COM 

3130 

PS I NT  96",  IL 

COM 

3140 

470 

PO  600  TCtl.TU 

COM 

3150 

IF  (APS  (COIFF  (L  EV.ITI)  .LT.l.E-2")  GO  TO  ECO 

COM 

3160 

MAXrIFOGM0>  «UF  V,  IT) 

COM 

3170 

PO  AFT  1:1, MAY 

COM 

3180 

J:MAXM|  -T 

COM 

3190 

M=ICI«  (LEV, IT  ,1) 

COM 

3200 

xa>=rTP(NH 

COM 

3210 

NII  =  IFX°(LFV,IT,I) 

COM 

3220 

22 


f  -r  — 


m*n 


4°  C 


se¬ 


gno 

90S 

91C 
92r 
93r 
94C 
950 
9sr 
9?  r 
98  0 
9Sr 


HiXF  CIi*MOJC(MI)  COM 

X?(J)=X(T>  COM 

COM 

TF  (  (MA*.  "'.I )  .09.  (IHtPM.ro  -1)1  GO  TO  490  COM 

PRII'T  940,  (HF XF< IP) ,HP ,IP=1 .MAX)  ,H  P  ,  H  P  ,H  ir  ,  (HFX  P  2  (  I P  )  ,H6  ,  I  P=  1  ,M  A  X  )  CO  M 
C  =  CriFFF(L‘VfTT)/?,  COM 

CALL  FR»rT  (C,IC,FT>  COM 

IF  (Nr.EO.F)  FPINT  c  1  0  ,  IC, MO, HIP, HP,  (HP,HP,IP=1,MAX),HB,MFL,  COM 

1  (HP, HP, TP^j,  MAX),  HR, HPp  COM 

IF  (  (ND.r.T.O)  . ANO.()O.LE.64>  1  PRINT  °?0,  TC  ,NP ,  HP ,  HI  F ,  M9  ,  (  NP,  HP ,  COM 
1 1°*l  ,MAX)  ,  1  ".FPL  ,  (HP, HP,  IP  =  1  ,MAX)  ,HP,HRP  COM 

TF  (N0.CT.F4)  PRINT  q»0,  C  ,H  a  ,HLP  ,H  p  ,  <MP  ,HB  ,  I°=  1  ,  MA  X  1  ,  H8  ,MFl  ,  COM 

1  (wp, hp, IP-i  ,  ha  x ) ,  hr  ,hr  p  COM 

P  =  INT  94-  ,  <  V(  IP)  ,HF,TP=  1  ,M»T)  ,H8,HP,  HP,  (v?  (IP)  ,h9,IF=1,M»X)  COM 

ro  TO  59"  COM 

P’INT  949,  (HfXP(IP),HP,I°=l ,  M  A  X  )  COM 

C3LI  FXS'T  (CpFFF  (MV, IT)  ,ir,ND)  COM 

IF  (NC.£".P)  FPINT  cio,  lC,MO,HB,H8,(HP,HP,IP*l,M»X)  COM 

TF  (  (  NO.FT  ,C  )  .  ANC.  (FO.  LE.64  )  )  PRINT  9  2  0  ,  IC  ,NC ,  H8 ,  HE  ,  t-8  ,  ( HP,  HO,  COM 

1 T°= 1  ,MAX )  COM 

IF  (Nr.GT.F4)  PRINT  930,  CO  r  Ff  (  LP  V,  IT  )  ,  H  P  ,  H  E  ,  H8  ,  (HP  ,  HE  ,  I  P=  1  ,  MA  X )  COM 
PRINT  940  ,  IT(IF)  ,HP,IP=l,MAx)  COM 

roNT  INO 7  COM 

f-0  TO  79  COM 

COM 

FCR"AT  <  T-»,7Ff  .1  ,n  (IX  ,»1,IZ  ,2V))  COM 

FORMAT  ()»M  .*••*•••*•  CMC  TIME  IS  NCN  s,F9.T,19H  EECONCS  •••♦•♦•COM 

1  •  •*  5X,»  1  9  /)  COM 

FCR-AT  ( 1  4  X  ,  If  ,f  0»1.)  COM 

FORMAT  (1  1  X,If  ,*■/*',  12,8041)  COM 

FORMAT  II 'y,F12.F,8i41)  COM 

FC»MAT  ( ?f  y ,8. 4 1 )  COM 

FORMAT  <1  up,->i.X“FO()ALG"/>  COM 

FCRM4T  (inr/lPO ,24X ,12/23* ,5H*IM  •/)  COM 

FORMAT  (“  f  POF  TCOIIPOSMAXILF.-W,  ITE  RM>  =",  3  16)  COM 

FORMAT  (-  fPOFTFOlITEPM2*",3I6)  COM 

FPPMAT  (“  APOF  TECITTERHL  (CEV>*",3I6>  COM 

FNO  COM 


3230 

3240 

3250 

326, 

3270 

3280 

3290 

33  00 
3311. 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 

34  30 
3440 
3450 
3460 
3470 

34  83 
3490 
3500 

35  1C 
3520 
3530 
3540 
3550 
3560 
35  7  0 
358C 
3590 
3600 
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SUPFOUf  TNC  r  A  N  C  FL  (1,1  TERM , I  TE,  ITD 

CAN 

10 

c 

TNI'  3U°R0UT  T,(S 

CAN 

14 

c 

A  IN  T  MIL  A  TFS  A  TFR«  WHOSE  COEFFICIENT  IS 

LESS  THAN  l.E-12. 

CA  N 

15 

CCM“CN  t"IF  <2,  199,11)  ,IEX°(2,199,li> 

, COEFF (2,199) , IF OSMAX (2,199) 

CAN 

20 

tf  (AncirnrFFll  ,TTE)).GT. 1.--1ZJ  RETUPN 

CAN 

30 

ITEFM=ITr®E  -1 

CAN 

40 

1TI  =  TTL-1 

CAN 

50 

TF  (ITt.LT  *ITr)  RETURN 

CAN 

SC 

P O  10  rTFTTf.TTL 

CAN 

70 

Na*  =  irrs«'ycL,iT*i) 

CAN 

60 

CORFF<l,TT)='*rFFF<L,IT  +  l) 

CAN 

90 

T  =  OSMA(  (L  ,  IT)  =  MAX 

CAN 

100 

no  in  rpFi,«sy 

CAN 

lit 

T0IR(L,It,IP)f1CIR(L,It  +  1,TP) 

CAN 

12  0 

It 

T2)(R(L,TT,TP)rIFXPIt,TT*l,IR) 

CAN 

130 

RETURN 

CAN 

140 

CNC 

CAN 

15  U" 

SUPPOUrTN-  FRACT  <C.,IC,,NO> 

C  THIS  «U°ROUTIMr  ATTrNPTS  TO  WRIT*  COEFFICIENTS  AS  INTEGERS*  OR 
C  FRACTION'S  TF'  T»FV  AFE  AT  LEAST  MULTIPLES  CF  1/64. 

NO=r 

TC=c*srsNt  i  ,f-cp * c> 
ip=vr 

CTFrOslPSCyr-FLC^TdC)  » 

IF  IDIF^O.  I  T.1  .E -OR)  RETURN 
rn  1"  1=1*7 

no=2**( 

riF=FLPTT(KP)«riFro 
TPiF=nrrn  *E-f  p 

IF  (OTF-cLfA T  (  I r I F>  .  LT.  1.5-.  6)  GO  TO  2C 

in  nuTiNU- 

PTCRM 

?o  ir=fo*roF  tsi'n(  ioif  ,ir> 

RETURN 

CNP 


FRA  10 
FRA  15 
FRA  16 
FPA  20 
FRA  30 
FRA  40 
FRA  50 
FRA  60 
FRA  70 
FRA  60 
FRA  °  0 
FRA  10  0 
FRA  11J 
FRA  120 
FRA  13  u 
FRA  140 
FRA  150 
FRA  160- 


Appendix  B 


EXAMPLES 


R1  COMMUTATORS  ARISING  FROM  DIAGON ALIZATION  OF  HAMILTONIAN 


[{Px.PzH.Py] 


Input  cards: 

2 

3  0.5  -0.5  XI  71  Y  ? 
3  0.5  -0.521  XI  Y  Z 

Result: 

112  2  11 

1/2  (  P  P  P  -  O  f  P  ) 

X  Z  Y  v  7  X 

11?  ’ll 

1/  2  (  P  P  P  -  P  F  °  ) 

Z  X  Y  Y  X  Z 

F0U»lS 


1 

♦  If*  * 

2  1  1  ’ 

1  (  p  o  «p  p  | 

X  Y  Y  X 

12  2  1 

-1  <  P  P  ♦  P  P  ) 
Y  Z  ’  Y 
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Input  cards: 


4  0.25  -0.’5  X  ?  v  ?  XI  xi 

4  0,?5  -11.25  X  ?  Y  2  XI  XI 

4  0.25  -0.25  Y  2  X  2  XI  XI 

4  0.25  -0.25  Y  ?  X  2  71  X  5 

Result: 


THE  TIME  TS  NOW  =  L.2F3  SECONCS  ***♦.*»»**  CYCLIC! 

7  2  11  It?? 

1/4  (  F  P  P  p  -  F  °  P  F  ) 

X  Y  X  7  7  X  Y  X 

2  2  11  112  2 

1/4  (FPFF  -  F  P  P  P  ) 

X  Y  2  X  X  7  Y  X 

2  2  11  112  2 

1/4(FCPP  .FPPP) 

Y  X  X  7  7  X  X  Y 

2  2  11  112? 

1/4  (  P  °  P  P  -  FPPP) 

Y  X  7  X  X  x  x  Y 

cOII4LS 


1 

♦  It*  * 

2  3  x  ? 

1  (  p  p  ♦  p  p  ) 

X  Y  Y  X 

4  1  14 

-1  (  F  P  +  F  o  ) 

X  Y  Y  X 

2  12  2  1  ? 

l(ppp  »pro) 

X  Y  7  7  Y  X 


7 

♦  IM  » 

2  1  1  2 

1  (  F  R  ♦  D  P  1 
X  Y  Y  X 

7 

1/2  p 

Y 

12  2  1 

-1  <  c  P  +  c  p  ) 

Y  7  x  y 
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R2  REDUCTION  OF  SINGLE  TERM 


P  P  P  P  P  P 

x  y  x  y  z  z 


Input  cards: 

1 

6  1.0  0.0X1  Y 


XI  f  1  Z  1  7  1 


Result: 


THE  TIME  IS  NCW  =  4. <.07  SECONPS  **********  CYCLICl 

111111 
1  P  P  P  P  p  P 
X  Y  X  Y  7  7 

EQUALS 

?  ?  7 
1  F  P  P 
X  Y  7 


1 

+TM  * 

1  1  3 
1  P  P  P 
X  Y  Z 


7 

♦  IM  * 

2  2 

1  F  P 
X  7 
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B.3  COMMUTATOR  ARISING  FROM  TRANSFORMATION  TO  MINIMAL  SET 


x  y  z 


z  y  x 


Input  Cards: 
4 


5 

0.25  -0.25  X 

1 

Y  1  7  1 

5 

0.25  -0.25  X 

1 

Y  1  7  1 

5 

0.25  -0.25  Z 

1 

VI  XI 

5 

0.25  -0.25  Z 

1 

Y  1  XI 

Result: 

1112 

2 

2  2  11 

1/ 

4  (  P  P  F  p  o 

- 

D  p  O  p  P 

X  Y  2  Y 

Y 

V  X  7  Y 

1112 

> 

2  2  11 

1/ 

4  (  F  P  P  p  ^ 

- 

p  p  p  o  p 

X  Y  Z  Y 

X 

Y  Y  7  Y 

1112 

2 

2  2  11 

1/ 

4  (  P  P  p  p  p 

- 

D  p  p  p  p 

7  Y  X  Y 

Y 

Y  X  X  Y 

1112 

2 

2  2  11 

1/ 

4  (  F  P  P  P  P 

- 

o  p  p  p  p 

7  Y  X  Y 

X 

X  Y  X  Y 

X  2 

X  2 
<  ? 
v  ? 


v  2 
x  2 

Y  2 
X  2 


50UAI.S 

1 

+  IM  * 

2  4  4  ’ 

-1  (  P  P  ♦  P  P  ) 
X  Y  Y  Y 

4  2  2  4 

1  (  P  P  +  P  P  ) 
X  Y  Y  Y 


+  IM 


2  P 

X 

2  2  2  ’ 

-5  <  P  P  ♦  P  P  ) 
X  7  7  X 

4 

-2  P 

Y 

2  2  2  2 

5  (  P  °  +  P  P  ) 

Y  7  Z  Y 
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c 

+X  M  * 


-8  c 

X 

2 

8  o 

Y 


R4  DIPOLE  MOMENT  COMMUTATOR 


[pr  •*' 

Input  cards: 

1 

2  1.0  -1.0  Y  2  U  1 


Result: 

2  1  1  *» 

1  <  r  p  -CD) 
YU  J  y 

fouais 


'  1 
♦  T  Y  * 

11  11 

1  (  P  P  ♦  P  o  ) 

W  Y  YU 


B.5  TRANSFORMING  INTENSITIES 

[i(PPP  +PPP)  IP2,  »  }] 
2  x  y  z  z  y  x  *  2  ^  x*  z  ' 

Input  cards: 


<• 


5 

0.25  -0.25 

X 

1 

Y 

1 

7 

1 

y 

? 

w 

1 

5 

0.25  -0.25 

X 

1 

Y 

1 

7 

1 

rt 

i 

X 

2 

5 

0.25  -0.  25 

7 

1 

Y 

1 

X 

1 

y 

2 

w 

1 

5 

0.25  -0.25 

Z 

1 

Y 

1 

X 

1 

w 

1 

X 

2 

29 


. 


— . . .  ■ 


■ 


IMF  TIME  IS  NCH  =  U.7J1  SFCONDS 


CYCLIC! 


1112  1 

12  111 

1/ 

4 

( 

P  P  P  P  3 

■  o  p  p  p  p  ) 

X  Y  7  x  W 

W  X  7  Y  X 

11112 

2  1111 

1/ 

4 

( 

F  P  P  P  P 

-  p  P  F  P  P  ) 

X  Y  Z  W  X 

X  W  7  Y  X 

1112  1 

12  111 

1/ 

4 

( 

F  F  F  P  n 

.  D  p  p  p  p  ) 

X  Y  X  X  A 

W  X  X  Y  7 

11112 

2  1111 

1/ 

4 

( 

P  P  F  P  a 

-  o  p  p  o  p  ) 

7  Y  X  W  X 

X  W  X  Y  7 

PQUflLS 


♦  IM 


12  2  •»  ?  1 

_1(opp  ♦  3  F  p  ) 

W  X  Y  v  x  W 

12  2  7  2  1 

1  <  P  °  P  f  3  F  =  ) 

W  Y  Z  7  X  W 

13  1  1  •»  1 


1/ 

2 

( 

PDF 

f  a  r  o  ) 

LX? 

XXL' 

1  2  1 

1  112  1 

1/ 

2 

( 

o  p  p 

p  *  r  o  p  p 

V  X  Y 

2  X  V  x  V 

7 

*TM  * 

ill 

1  1  1 

- 

3 

( 

F  F  F 

»  F  F  P  > 

U  X  7 

Z  X  U 

111 

1  1  1 

2 

( 

F  c  F 

t-  P  F  o  ) 

V  Y  Z 

T  Y  V 

1  2 

’  1 

3/ 

2 

( 

P  p  ♦ 

P  P  ) 

M  V 

Y  W 

1  2 

2  1 

3/ 

2 

( 

P  o  + 

P  P  ) 

W  7 

x  w 

